Critical Scaling of Extended Power Law I-V Isotherms near Vortex Glass Transition 
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In view of the question about the vortex glass theory of the freezing of disordered vortex matter 
raised by recent experimental observations we reinvestigate the critical scaling of high T c super- 
conductors. We find that dc current- volt age characteristic of mixed state superconductors has the 
general form of extended power law which is based on the Ginzburg-Landau (GL) functional in the 
similar way as the vortex glass theory. Isotherms simulated from this nonlinear equation fit the 
experimental I-V data of Strachan et al.[Phys. Rev. Lett. 87, 067007 (2001)]. The puzzling ques- 
tion of the derivative plot for the I-V curves and the controversy surrounding the values of critical 
exponents are also discussed. 

PACS numbers: 74.25. Sv,74.25.Fy, 74.25. Ha,74.20.De 



The hypothesis of a vortex glass in disordered high 
temperature superconductors |l| has spurred much re- 
search and many discussions during more than one 
decade and continued to be a very controversial issue 0, 
0. Though the striking dc current- voltage (I-V) data 
collapse shown by Koch et al. 4], and found by other 
groups on other samples, strongly supports the consen- 
sus that a glass transition occurs in high-temperature 
superconductors(HTS) 5], this widely accepted evidence 
has met serious doubt recently. Strachan et al. show wide 
range accurate isothermal I-V measurements over 5 or 6 
decades on a high quality 2200 A thick YBa 2 Cu 3 7 -x 
film with T c w 91. 5 K and transition width about 0.5K|3j 
and find although the I-V isotherms measured in a mag- 
netic field can be collapsed onto scalling functions pro- 
posed by Fisher et al. [lj as is widely reported in the litera- 
ture, these excellent data collapse can also be achieved for 
a wide range of exponents and glass temperatures T g as 
demonstrated in their Fig. 2 (a), 2(b) and 2(c) 3]. Since the 
critical temperature T g cannot be determined uniquely, 
they argued that necessary evidence for a vortex-glass 
transition to zero resistance is still not seen in their I-V 
experimental data, as well as in all the ones they know of 
in the literature which scale. Therefore, the correctness 
of the vortex glass picture has to be reinvestigated. In 
this work, we try to settle this remarkable controversy 
by introducing a general extended power law form of the 
current- volt age isotherms. 

In analogue with the vortex glass theory Q, our ex- 
tended power law equation is also based on the GL free 
energy functional of a system of N flux lines with a field 
H along the z direction in a sample length L which can 
be described with the free energy represented by the tra- 
jectories {rj (z)} of these flux lines and the pinning po- 
tential Vp (r) arising from inhomogeneities and defects 
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Thermally activated flux motion can be considered as 
the sequence of thermally activated jumps of the vor- 
tex segments or vortex bundles between the metastable 
states generated by disorder. Every elementary jump is 
viewed as the nuclear ation of a vortex loop, and the mean 
velocity of the vortex system is determined by the nucle- 
aration rate Q, Q 



v oc exp (SF/kT) 
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Here SF is the free energy for the formation of the critical 
size loop or nucleus which can be found by means of 
the standard variational procedure from the free energy 
functional due to the in-plane displacement u (z) of the 
moving vortex during loop formation 
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where /l = J x e^/c is the Lorentz force due to applied 
current J and f v is the viscous drag force on vortex, 
frj = -rjv VO rtex, with v vor tex = d~u / dt and viscous drag 
coefficient 77 ~ (<&oB C 2) / (p n c 2 ) as estimated by Bardeen 
and Stephen p. 

There are two essential differences between our equa- 
tions (1) and (3) and the corresponding ones used in pre- 
vious vortex glass theories: 

(i)We considered both the Lorentz force / l and the 

viscous drag force / v on the vortex loop during its for- 
mation while in present vortex glass theories only bare 
Lorentz force / l has been taken into account 0,0- 
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(ii)The effect of the real sample size L on the current- 
voltage response is explicitly estimated in our case. The 
influence of the finite size effects on the I-V isotherms is 
recently revealed even in thickest films at zero magnetic 
fieldQ. 

These two factors result in a new form for the I-V 
isotherms of superconducting vortex matter. 

In stead of the barrier energy found by the vortex glass 
model IK as 



SF = U(J) « U c (^fr 
now we have the barrier energy 

SF = U(J P ) « u c (^r 
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which implies a current- voltage characteristic of the form 



E(J)=p f Jexp[-^(fr] 



(7) 



where U c is a temperature- and field-dependent charac- 
teristic pinning energy related to the stiffness coefficient 
and J c is a characteristic current density related to U c ,/i 
is an numerical exponent. 

Considering the real size effect Eq.(7) leads to a general 
normalized form of the current- voltage characteristic 
in the form 



y = xexp [-7 (1 + y - xf] 
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with Jl the transport current density corresponding to 
the case where the critical size of loop formation is equal 
to the real sample size L. The pf in Eq.(7) is the flux 
flow resistance of a pinning free mixed state as derived 
by Bardeen and Stephen. Its relation with the normal 
resistance is of the form 
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In an earlier work [TT[, equation(8) has also been shown 
in connection with the Anderson-Kim model. Com- 
pared with the widely used power law nonlinear response 
E oc J n (see Ref.5), equation(8) can also be nearly equiv- 
alently expressed in an extended power law form 
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with n a numerical factor depending on the parameters 
7 and p in Eq.(8), which can be determined from solving 
the equation 
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FIG. 1: Comparison of solutions of Eq. (8) (solid symbols) and 
Eq.( 10) (curves). 



where the subscript i denote the corresponding values at 
the inflection point of Iny ~ Inx curve for Eq.(8) with 
the maximum derivative S max so 
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In Fig.l we show the numerical solutions of Eq.(8) and 
Eq.(10) for comparison. 

The complex experimental current-voltage data are 
usually described phenomenologically with a formQ 



E(J) = J Pf (T,B,J)=Jp f e 
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Equation(8) implies that effective barrier can be explic- 
itly expressed as 

U eff (T,B,J) = U c (B,T)F[J/J c (B,T)} (13) 

withJ c (£,T) = J L (B,T). 

Incorporating it into the commonly observed scaling 
behavior of magnetic hysteresis M(H) in superconduc- 
tors, it can be shown that U C (B,T) and J C (B,T) in 
Eq.(13) must take the following forms |5| 

U C (B,T) = ^(T)B n ex [T*(B)-T] 5 B n 
J C (B,T) = X(T)B m ex [T*(B) -T] a B m (14) 

with T*(B) being the irreversibility temperature where 
tilt modules and stiffness of vortex matter vanish. Con- 
sidering pf ex T, one finds from Eq.(8) the following tem- 
perature dependent relations for 7 and I — V isotherms 
of optimally doped YBCO samples in a given magnetic 
field as 

7(T) = 7o (T*-T)V^T 

IocxJ L (T) ocx(f -T) a 

V <xyJ L (T) Pf <xy(T* -T) a T (15) 
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FIG. 2: Comparison of simulated I-V isotherms 
from Eqs.(8) and (15) with the experimental data of 
Ref.3[Phys.Rev.Lett.87, 067007(2001)]. 

(a) I-V isotherms simulated from the extended power law 
equations (8) and (15) with assumed T*(B = 4T) = 88i^T, 
and a = 3.0,6 = 1.5, p = 1.75 and 7o = 3.21 x 10 21 . 

(b) Comparison of our simulated I-V isotherms (lines) with 
the experimental data of Ref.3(solid symbols). 



Combination of equations (8) and (15) enables us to sim- 
ulate the experimentally measured I — V isotherms of 
Strachan et al. 3] as shown in Fig. 2, where the param- 
eters used in (8) and (15) are a = 3.0,(5 = 0.5, p = 
1.75,7o = 3.21 x 10 21 and the irreversibility tempera- 
ture T*(B = 4T) is assumed SSK similar to the Fig.l 
in Ref.[3]. On the basis of nearly 100 simulated I — V 
isotherms from 6 8 if to 8 8 if with temperature intervals of 
0.2if , we obtain three different scaling data collapse anal- 
yses with assumed critical temperature T g of 81if ,75if 
and 70if respectively. In Fig. 3 we show the comparison 
of our result with the data collapses of measured I — V 
curves of Strachan et al. in Ref.[3] and find fair agree- 
ment. 

In Fig. 3 of the Ref.[3] another serious puzzle is men- 
tioned with plotting the derivatives of the InV vs Inl for 
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FIG. 3: •: Collapse of the simulated I-V isotherms in Fig. 2(a) 
using various critical parameters compared with 
o: the experimental I-V isotherms for a 2200 A YBCO film in 
4T measured by Strachan et al.[3]. The conventional analysis 
is shown in (a). 



experimental isotherms following Repaci et al. 12]. All 
isotherms seem to have a maximum slope at about some 
current with Ohmic tails developing to the left of the 
peaks. This behavior is inconsistent with the scaling 
function of vortex-glass transition 4] which should pre- 
dict a horizontal line in the [dlnV/dlnI]T ~ Inl plot at 
T = T g . Nevertheless, our extended power law I — V 
equations may give it a reasonable explanation. 

In Fig. 4 we plot the slope of experimental I — V 
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FIG. 4: Empirical scaling of the derivatives of the experimen- 
tal InV ~ Zn/ isotherms. A:data from Ref.12 B: data from 
Ref.3. The solid curve is the results simulated from Eq.(17). 



isotherms of Ref.[3] and [12] in an empirical scaling form 

= *(//*) (16) 

and find that these experimental data are collapsed, 
where the subscript i denotes the value at inflection 
point of InV ~ Inl curve, S = [dlnV / dlnI]T and n de- 
notes Smax at inflection. This empirical scaling function 
can be derived from extended power law nonlinear I-V 
isotherms. Starting from our equation (8) one finds ana- 
lytic equation 

S-l = c (x-y) In (x/y) 
n — 1 1 + y — x + yp\n (x/y) 

where 

c = l + Vi-Xj+yiphi (xj/yj) 
(xi - yi) In (xi/yi) 



shown by the solid curve in Fig. 4. 

The agreement of the simulated I-V isotherms derived 
from Eqs.(8) and (15) with the wide-range experimental 
data of Ref.[3] as shown in our Fig. 2, 3 and 4 demon- 
strates that extended power law is the general form of 
I-V isotherms. Though, in principle, each I-V curve can 
have both positive and negative curvature parts joined 
at the inflection point, the actually used current- voltage 
window in measurements limits the observable part of 
each I-V isotherms and results in the previous reasonable 
scaling analysis Q, Q . Using wider I-V windows of mea- 
surements, Strachan et al. find that non-Ohmic-power- 
law-like behavior can be fit to all T < 81 if isotherms 
over at least four decades of voltage data and thus in- 
sightfully reveal that the conventional I-V scaling can- 
not determine a unique Tg with corresponding critical 
exponents. The successful reproduce of their experimen- 
tal data collapses with our simulated I-V isotherms from 
equations (8) and (15) with the uniquely assumed param- 
eters the irreversibility temperature T*(B = AT) = SSK 
and a 3.0, S = 0.5, p = 1.75, 70 3.21 x 10 21 shows 
that extended power law is helpful in settling the contro- 
versy surrounding critical phenomena in high tempera- 
ture superconductors. 

In conclusion, we show that dc current- voltage charac- 
teristic of mixed state superconductors has the general 
form of extended power law. Isotherms simulated from 
this nonlinear equation fit the recent experimental I-V 
data of Strachan et al. This equation also gives a reason- 
able explanation to the critical scaling relation recently 
observed in the derivatives of the InV vs Inl. 
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